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ABSTRACT 

Interstellar matter that is accreted onto isolated magnetic neutron stars in the Galaxy 
(~ 10^ by number) is accelerated and reflected back by MHD shocks, which envelope the 
stars. The integrated power in the Galaxy Lcr,ns is ~ lO^'^erg s~^, the energy distribution 
is a power law of spectral index > 2, and the particle energy can be raised to 10^ GeV, 
consistent with the power and spectrum of primary cosmic rays in the Galaxy. 

The major contribution for Lcr,ns comes from a minority of ~ 10^ isolated neutron 
stars which are located within dense clouds. Sources in these clouds, that are generally 
spread within the Galactic disk, can explain the concentration of high-energy cosmic rays 
in the Galactic plane, as deduced from pion decay spectra in gamma-ray observations. 
The soft X-ray luminosity from these neutron stars is consistent with the Galactic X- 
ray background. The accretion may be associated with ion-neutral bias, that is further 
enhanced by ion conflnement in frozen-in magnetic flelds, which can raise the relative 
abundance of flrst ionization potential (FIP) elements in the cosmic rays. 

Key words: cosmic rays, acceleration of particles, shock waves, stars: neutron. X-rays: 
ISM 
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1 INTRODUCTION 



The origin of cosmic rays has long been a central question in astrophysics (e.g. Cesarsky 
1987, hereafter C87, Blandford & Eichler 1987 [BE87], Axford, 1994). The angular dis- 
tribution of cosmic rays is highly isotropic. The isotropy accuracy is ~ 0.1% at energies 
< 10^ GeV and ~ 1% at 10® GeV, but it is limited to a few tens of per cent in the ultra- 
high-energy range > 10^ GeV. The origin of cosmic rays, apart of the ultra-high-energy 
class, is believed to be Galactic, associated with continuous leakage from the Galactic 
magnetic field confinement. The confinement time is observationally constrained by the 
decay time of unstable nuclei, e.g. ^°Be {td — 2.2 x 10^ yr), or ^*^A1 {ta — 0.85 x 10^ yr), 
to be ~ 8 X 10^ yr (for homogeneous ISM). 

Three major requirements of candidates for cosmic ray sources are that they produce 
the observed cosmic ray power, spectra, and elemental abundance. The total power of the 
Galactic component is estimated to be Lcr 2 - 3 X 10^° erg s"^ (BE87). The spectra, 
which are noisy and depend on the specific ion, are approximately a power low oc E~'^, 
where a < 2.5 at energies below ~ 10^ GeV, ~ 2.7 up to 10^-^ GeV and ~ 3 at > lO*^-^ 
GeV. The elemental composition (Meyer 1985a,b, Silberberg et al. 1987) at energies up 
to 10^ GeV clearly deviates from the Local Galactic abundance, but fairly agrees with the 
Coronal composition. Specifically, there is a relative enhancement in cosmic rays of the 
low FIP (< 10 eV) elements. 

Supernova remnants (SNRs) that heat and accelerate the surrounding material via 
shock wave acceleration are widely accepted as the origin of Galactic cosmic rays (Bell 
1978, Blandford & Ostriker 1980, BE87). The cosmic ray power is only a few per cent 
of the total power associated with SNR mass motions, therefore, even an inefficient shock 
acceleration mechanism is applicable. Such a mechanism also explains the spectral shape. 
The maximal energy gain per particle E^ax, however, is limited by the shock lifetime, 
estimated to be ~ 10^S/(10~^G)GeV in optimal conditions {Emax may be increased to 
5 X 10^ GeV if the SNR is propagating into a stellar wind cavity [Volk & Biermann 1988]). 
Moreover, the FIP bias is unexpected, since the SNR is accelerating material which is 
believed to have the Local Galactic abundance (this problem is relaxed in models of SNR 
in a stellar wind cavity, if the presupernova object is a red or a blue giant, or a Wolf-Rayet 
star [Silberberg et al. 1990]). In competing 'stellar' models sub-relativistic particles are 
injected from F and G stars (Meyer 1985a,b), which explain the cosmic ray composition, 
but require a secondary stage of acceleration, e.g. by passing a supernova shock wave. An 
injection mechanism for ultra-high-energy cosmic rays was proposed by Kazanas & Ellison 
(1986), based on accreting binary systems such as massive X-ray binaries, where a strong 
stellar wind of the companion star is accreted onto the compact object, a neutron star 
(NS) or a black hole. 

Cowsik & Lee (1982) have shown that accretion fiows with shocks onto isolated NSs 
can accelerate cosmic rays very efficiently. They consider old NSs, assuming that their 
surface magnetic fields are weak and do not affect the hydrodynamics. The fiow stagnation 
imposed by the stellar surface results in a standing spherical shock wave which envelops 
the star very close to its surface. Potential energy, associated with interstellar matter 
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(ISM) accreted onto the star, is continuously channeled to acceleration of high energy- 
particles. The cosmic ray power is ~ 10^^ erg s~^ per NS, assuming that its velocity is 
smaller than the ISM sound speed, and the ISM density is large ~ 10^ cm~^. For such a 
luminosity a few times 10*^ NSs are sufficient to explain the Galactic cosmic rays. Such 
a high power, ~ 10~^ of the Eddington luminosity, may be, however, associated with a 
temperature increase above the pair production threshold, as was recently shown by TuroUa 
et al. (1994). Pair annihilation and gamma ray emission can significantly cool the particles 
and (electron) synchrotron would become very efficient, whether the magnetic fields are 
originated in the NS or frozen within the accreted ISM. Also, if the peculiar velocity of 
an old NS is larger than the ISM sound speed the available power is significantly smaller 
than 10^^ erg s~^. More importantly, various observations suggest that the magnetic field 
of an old NS remains very strong (e.g. Srinivasan 1989, Phinney & Kulkarni, 1994), and in 
that case the magnetic pressure can halt the fiow at the Alfven surface far from the stellar 
surface. 

In this paper I focus on Galactic cosmic rays that originate from accreting magnetic 
isolated NSs. Because of the halting of the accretion, a standing shock is created, within 
which shock acceleration occurs (Shemi 1995). The available power, compared with the 
free fall luminosity, is significantly reduced, to a value orders of magnitude below the 
Eddington luminosity. Other mechanisms that might regulate the accretion, such as the 
"propeller effect" (lUarionov & Sunyaev 1975) or a relativistic wind induced by the NS, 
which will not be discussed here. Accretion through the polar caps, which, unlike disk- 
accretion, probably is inefficient in spherical inflows, is ignored, and instabilities of the 
magnetosphere are not considered. The accreted matter is never heated to ultrarelativistic 
temperature, and pair production therefore does not take place. 

The total population of Galactic isolated NSs is estimated to be A?" ~ 10^ (Blaes & 
Madau 1993, hereafter BM93), but the major contribution to cosmic rays comes from a 
small subgroup (~ 1%) containing the slow NSs which are located within dense clouds. 
In the following paper I first describe the physics of accretion, shock waves and particle 
acceleration around an isolated magnetic NS, and then discuss the applicability of these 
objects as sources of cosmic rays. 



2 ACCRETION ONTO A MAGNETIZED NEUTRON STAR 

An isolated NS that is propagating with velocity v^q = v/(40 km s"-*^) accretes ISM from 
a cylinder of a Bondi-Hoyle radius 

2GM 

TBH = -^-^ ^ 2.3 X lO^^v^o'^ cm . (1) 

Here Cg is the sound speed and M is the NS mass, where throughout this paper we 
use M = 1.4M0. I will focus on dense ISM clouds (n = 10^^77.4 cm~^), where the 
accretion rate would be large. Extreme UV and soft X-ray radiation emitted by the 
accreted plasma is sufficiently intense to ionize the matter within a Stromgren sphere 
~ 3 X 10^^(S'/10^^)n4 ^^^(cKs/lO"-^^)"-^/^ cm, where S is the rate of ionizing photons. 
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in photons s"-*^, and cub is the recombination coefficient to excited states of hydrogen, in 
cm^ s~^. In the ionized gas macroscopic magnetic fields and collective plasma phenomena 
couple the particles, and the accretion become hydrodynamical, rather than coUisionless. 
Removal of angular momentum becomes efficient and raises the accretion efficiency, by a 
factor as large as {c/csY ~ 10^, compared with that of a neutral gas (e.g. Shapiro & 
Teukolsky 1983). 

Far from the central object (r* ^ r < tbh, where r* is the NS radius, 10^ cm in this 
paper) irreversible processes of energy dissipation are negligible, and the pressure obeys 
the adiabatic relation 

p=po(-)^ (2) 

where F is the specific heat ratio. The expected accretion rate is 

M = AtttIh^sPv ~ 4.3 X 10^V^^4o^g ^"^ , (3) 

where A., is a numerical factor of order of unity. The accretion rate is much smaller than 
the Eddington rate MecI = LEdT*/ {GM) ~ 10^^ g s~^, therefore radiation pressure is 
never strong enough to halt the accretion. 

The continuous removal of angular momentum causes streamlines to bend and con- 
verge toward the NS. To follow the flow dynamics we neglect the influence of magnetic 
fields and back streams from the central object or from the shock, and assume spherical 
symmetry. The dynamics are governed by the continuity equation (3) and by the Bernoulli 
equation 

\ = . (4) 

2 F- 1 r ^ ' 

The sound speed [c^ = {^ikT jmY^I'^ = 0.57 (T/50 K°)V2 km s-^] of the dense and cold 
ISM cloud is low, therefore, even for small NS velocity, the flow at r = vbh is supersonic, 

- —1/2 

with a rather large Mach number M = 70 V4QTr^Q . Unless nonthermal acceleration 
occurs, the energy gain per particle does not exceed [1 — ("^//'/c)^]"^/^ — 1 = 0.34 of its rest 
mass, maintaining the plasma subrelativistic even when approaching the NS, with typical 
blackbody temperature ~ 10^ — 10^ K. Substituting the nonrelativistic value F = 5/3 and 
the limit Cs <^ v, equations 2-4 give the scaling behavior of v, Cs, p, and the gas ram 
pressure pv^: 

V oc r~^^'^ ; Cs oc r~^^'^ ; p oc r~^^'^-, pv'^ oc r~^/^ . (5) 

Note that the Mach number remains constant. The ionized gas has high conductivity 
and the magnetic Reynolds number is large, so that the residual magnetic fields in the 
ISM are frozen within the infiow. Field lines are swept along the streamlines and the 
global field pattern becomes radial [one may consider a large - scale field which initially is 
homogeneous {Bn^, By, B^) = (-Bo,0,0), in where an accreting object is placed. The lines 
of the frozen-in field would be bent towards the central object, along the y and z axis, 
resulting in a configuration of a nozzle. The nozzle configuration is narrow in its central 
part and globally the field lines surrounding the central object become radial]. The field 
strength increases oc and the magnetic pressure would therefore increase as r~^, faster 
than the ram pressure and the gas pressure (F~-^pCg oc r~^/^). Such a field enhancement 
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would, however, be regulated when small scale magnetic reconnections become sufficiently 
fast. The field strength is therefore bounded by the cquipartition limit, that is achieved 
when the reconnection time scale, t/va, becomes comparable to the compression time 
scale, r/vff^ i?^„g/(47r)~^ < pv^ . This occurs at the cquipartition radius 

= ^ ^ 10«(3|i^)"\r/=.-°'' cm (6) 

where Bjsm is the magnetic field strength at tbh- 

Magnetic fields of the strength expected in old NS, 10^ < -Bsurface ~ lO^^G (e.g. 
Srinivasan 1989, Phinney & Kulkarni 1994), would cause stagnation of the flow at the 
Alfven surface, where the magnetic pressure balances the flow ram pressure. The flow is 
decelerated and compressed, and standing shock waves occur ahead of the Alfven surface. 

Using the dipole approximation B oc r~^, the Alfven radius is given by 

Tm = lQ'B%\-'/'v%' cm . (7) 

Here Biq is the NS surface magnetic field in units of 10"^° G. 

If the light corotation radius r/c = cP/^tt = 4.8 x lO^PiS (Pi = P/1 s) is larger then 
Tm a pulsar wind beyond the light cylinder can halt the inflow. However, such a mechanism 
would be signiflcant if the NS rotation period is shorter than 

= 2 X lO-^sJ/'n^:'/'^^/" s , (8) 

while old NSs, unless recycled, are generally slower than ~ 1 s. Beyond ric the magnetic 
field decreases as r~^, therefore in very fast rotators (P < Pic,m) the Alfven surface would 
generally be pushed beyond tbh- If > tbh-, the accretion is prevented in the first 
place. The motion of the NS through the ISM is then associated with a stationary bow 
shock, distant 

Tism = 3.2 X 10^^Pio^;4o' Pf^ cm (9) 

ahead of the star. The accretion flow will overcome the centrifugal barrier (where the 
free fall velocity become equal to the Keplerian velocity), at the centrifugal radius Vcent = 
[GM/Vt^y/^ = 1.7 X IQ^pI'^ cm, if the star period P exceeds 2TTrU'^{GM)-^/'^ . With 
equation 7 this condition reads 

P > Pcent,m - QAQB%\f"v^'^ S . (10) 

For NSs in dense ISM clouds we can reasonably assume P > Pcent, m, and therefore 
consider a shock wave at Vm ~ 10 V r 8 cm. 

3 PARTICLE ACCELERATION 

The topology and physics of the shock ahead of the Alfven surface is complex and not fully 
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resolved yet. Shock topology would differ from a spherical shape, presumably is conical 
with an opening angle ~ and a bow front, and the obliquity varies from parallel 

to perpendicular. A large electron and ion concentration is built up in the foreshock 
regions. Irregularities in the shock topology and temporal variability are also expected, in 
the same manner as in simulations of wind accretion onto an unmagnetized massive object 
(Ruffert, 1994 and references therein). Backward streams and perturbations influence the 
topology, sometimes depending on the Mach number and the r_B///r* ratio, making the 
flow configuration unstable. 

The jump conditions for a planar-discontinuity, infinitesimally narrow shock imply 
a compression ratio ( = P2/P1 (here '1' denotes the upstream region and '2' stands for 
downstream) that approaches the limiting value ( ^ (F + l)/(r — 1) = 4 for F = 5/3 
(e.g. Draine & Mckee 1993). The parallel component of the magnetic field is conserved 
(i?i 11 = -82,11)7 while the perpendicular component downstream is enhanced (-82,^ = C-8i,± 
for a flow with no transverse velocity component on both sides). The true strength of the 
magnetic fields upstream is reduced if the NS dipole fields are screened by electric field 
[v X B) in the shock region. On the other hand, the enhancement of frozen magnetic fields 
in the inflow would increase the upstream fleld. 

In the following I will neglect the complexity in the shock conditions, and assume that 
known acceleration mechanisms in planar shocks are still applicable, to assess the cosmic 
ray luminosity, the maximal energy per particle and the energy distribution. 

Some energy loss mechanisms take place at the shock region, among them radiation, 
dissipation through MHD waves and escape of neutral particles. Deflning fcr the fraction 
of the free fall luminosity {Iff — MGMr~^) that is left in cosmic rays after such losses, we 
obtain the power per an accreting NS 

lcr,ns = fcrM^ = 8 X 10^ VcrA«i?ro'^'n^/ Vo'"^" erg s-1 . (11) 

Particles would be accelerated through multiple shock crossing, presumably via a dif- 
fussive shock first-order Fermi mechanism (Drury 1983, BE87). Shock drift acceleration, 
moreover, can take place in the quasi perpendicular segments of the shock, where an elec- 
tric field V xB exists. The nuclear mean free path Icoi = ^pp/nm,p (where Xpp ~ 60g cm~^) 
at the shock region is ~ 3 x lO^^n^^ cm, much larger than the shock scale r^. However, 
the proton Larmor radius Ib = ^mpC^/{eB) will be smaller than unless particles are 
ultra relativistic (7 ^ 10^). Particle - wave interactions will therefore dominate the inelas- 
tic (particle - particle) collisions, making the shock coUisionless. Particles are scattered 
by Alfven waves, small-scale magnetic irregularities and turbulences, or other collective 
phenomena. 

If the upstream magnetic pressure Pmag is well below the gas pressure Pgas, the MHD 
perturbations are basically Alfven waves that travel upstream with velocity va < vi. In 
the downstream region the plasma experiences strong magnetic turbulences. Particles 
are reflected by these turbulences backward to the upstream region, then reflected again 
downstream by the Alfven waves (or also by turbulences, if Pmag ~ Pgas upstream. Note 
that even in that case we expect a large pressure gradient which would force a shock). 
For a particle with velocity w the momentum gain per single shock crossing is p — > p x 
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4[v2 — vi]/3w. After a sequence of ~ c/w random crossings the momentum is raised to the 
rclativistic regime, and, since w ':$> va, the Alfven waves upstream are seen as a quasi-static 
barrier. 

When the acceleration is sufficiently rapid, and not affected by energy losses, the 
particle energy is raised until the particles diffuse from the shock. The conical shape of the 
shock allows downstream accelerated particles to escape. The diffusion coefficient may be 
assumed to have the general form k = Ato/3, where A is the particle mean free path, usually 
between one to a few Larmor radii. Particles remain confined within the shock as long as 
the upstream diffusion length scale k/vi is larger than the shock radius of curvature ~ rm- 
The maximal energy per particle can be estimated by using the Bohm approximation 
(A = Ib) and requiring (Ibw) /{3vi) < r^, that gives (with w = c; vi = {2GMy/^rm^^'^) 

Emax = Imaxmc" ~ 5.4 X lOV-^/'Sio GeV . (12) 

To obtain a value of Emax comparable to the 'knee' value < 10^ GeV, at which the cosmic 
ray spectrum bends, one needs rm,8 ~ 0.3. Since the Alfven radius scales 
Emax scales oc S""^/^, and such large values of Emax does not have to be associated with 
a large accretion rate, if the stellar magnetic field is small Biq <C 1. 

Energy loss by synchrotron emission also can, in principle, regulate the particle ac- 
celeration. The maximal energy in that case is determined by the competition between 
these processes. The time scale for diffussive acceleration tacc = p/p ~ ^/"i^i [more rig- 
orously, taking care for the time spent by a particle in each side of the shock, one can 
obtain tacc = 3 x (t^i — V2)~^{ni/vi -\- 1^2/^2), e.g. Drury 1983] while the ion cooling 
time scale tsynch,ion ~ 6 x 10^^A^Z~'^B~'^'y~^ s (here B in Gauss). From the requirement 

tsynch,ion > tacc WC obtain the Constraint ^max.ion ~ 

2.6 X 108yl3/2z-2s-'/V^,8 (notc 
that in this case ^ max, ion only weakly depends on B, namely ^max B^/^'^). 

The synchrotron cooling time of electrons is, however, shorter then that of protons, by 

a factor (me/mp)^, implying (with -fmax,e « im^/mpY/'^) ^imax^e ~ 3.3 x IQ^ B^Q^'^rm,8 S- 
Synchrotron cooling although negligible for ions below ~ 10^ GeV, becomes acute for 
electrons. Moreover, electrons acceleration is not even started unless additional 'injection' 
mechanism operates. MHD waves upstream do not resonate with electrons below ^min,e — 
1 = mp/rrie x {va/cY (for {c/vaY ^ ij^p/'n^eY^ hence only the relativistic electrons can 
be confined in the shock. The mechanism discussed here therefore may not be applicable 
for cosmic ray electrons. 

The initial energy distribution of particles upstream is roughly a delta function at 
GM/rm ~ (0.064c)2r-V with some thermal broadening {kT ~ <C v1). By the stochas- 
tic nature of the acceleration, and since the total number of scatterings per particle is 
sufficiently large, the distribution is spread out with no preferred scale. From the linear 
theory of acceleration in subrelativistic shocks we expect a general form 

n{E) oc . (13) 

The spectral index, given by 

«=^~2, (14) 
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is independent of the magnetic fields, and is found to be typical in other cases of shocks. 
Cowsik & Lee (1982) show that a power law typical to cosmic rays can also be obtained 
in shock waves which have a spherical geometry. Ellison & Eichler (1985) also show that 
q; > 2 in non-linear shock models that account for the cosmic ray feedback. Even when 
the Alfven Mach number gives rise to strong compression, C ^ 4, the slope is almost 
unaffected. 



4 DISCUSSION 

The Galactic isolated NS population (BM93) is estimated from the pulsar birthrate, as- 
suming a steady state, tohe N > 10^, but nucleosynthesis constraints on Galactic chemical 
evolution require a total number as large as A?" ~ 10^ (Arnett, Schramm, & Truran 1989). 
A fraction of the potential energy associated with matter accreted onto these NSs is contin- 
uously channeled to accelerate particles to very high energies. Although the detailed shock 
structure and acceleration process are yet unclear, the energy budget and the knowledge 
of acceleration in astrophysical MHD shock waves suggest that accreting isolated NSs are 
significant cosmic ray sources. The cosmic ray luminosity Lcr,ns is the integral over the 
cosmic ray emission lcr,ns from all the Galactic isolated NSs, and other Galactic compact 
accretors like isolated white dwarfs, black holes, and massive binaries, that should be taken 
into account. 

We briefly discuss the parameters that affect the accretion rate and thereby the cosmic 
ray luminosity, namely the ISM density distribution, the NS velocity distribution, their 
surface magnetic fields and periods, and the physics behind the efficiency factor f^r- 

The filling factor of the ISM dense phase is estimated to be ~ 1% in the disk, and 
it rises to ~ 10% in an extended region in the Galactic center. From the estimated NS 
density in clouds 5 x lQ-'^{nns/vc~^){Maoud/MQf/'^, BM93 deduced ~ 10^ isolated NSs 
in the 19 clouds closest (< 800 pc) to the Solar System. 

The rms value of their velocity distribution is calculated to be 98 km s~^, where a 
substantial fraction of them, fs ~ 0.25, are slow [v < 40km s~^), and fa ^ 0.2 of them are 
located within dense ISM clouds (n ~ 10^ — lO^cm"^). The upper density value represents 
the cores of giant molecular clouds. BM93 show a disk velocity distribution with 6% of 
NSs having v < 20 km s-^ 22% below 40 km s-\ and 50% below 72 km s-^ It turns 
out that the major contribution to Lcr,ns would come from a minority of 10^ isolated NSs 
which are slow and located within dense clouds. 

The strengths of magnetic fields of old NS are controversial. Field decay is predicted 
in many NS models, which have also attempted to explain why such decay slows down 
or stops. However, there are arguments against a decay mechanism, in which case the 
diversity in pulsar magnetic fields is the result of the conditions during their creation (see 
e.g. Srinivasan, 1989 Phinney & Kulkarni 1994). 

The fraction fcr of the input power that is channeled to cosmic rays depends on the 
efficiency of competing energy loss mechanisms. Synchrotron radiation can maintain the 
electron temperature low but only loosely affect the ions. Free - free emission [e// x 
~ lO^^n^T/lO^ K)-*^/^ erg s"-*^] and Compton losses are also small. Escape of neutrons. 
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which are produced through pp —>■ nX reactions or photodissociation of ^He nuclei, can 
significantly ( < 25%) afi'cct the ultra-high-energy part of the spectrum (Kazanas & Ellison 
1986). Energy is also dissipated through MHD waves. For the sake of evaluating Lcr,ns we 
will use fcr ~ 0.25 and focus on the slow stars in dense clouds. Using the canonical values 
discussed above we have 

L,rns = 1.2 X 104O(J?^)(^)(A)(i£!l)(:^) erg s'^ (15) 
' UO^^^O. 25^^0.2^^0.25^40^3^ & > v J 

which is consistent with the power estimated for Galactic cosmic rays. 

The maximal energy gain per particle (Equation 12) is sufficient to explain the major- 
ity of cosmic ray particles, but not the very-high-energy ones. With a large accretion rate 
Emax would greatly be raised, to values similar to those proposed by Kazanas &; Ellison 
(1986) for binary systems. Moreover, coherent acceleration through magnetic field anni- 
hilation (frozcn-in fields -|- NS fields) as is suggested here, as well as pulsar acceleration, 
as proposed elsewhere by many authors (e.g. Srinivasan 1989, Arons & Tavani 1994), are 
likely to take place. 

Of particular interest is the subgroup of very slow {v < Cs) isolated NSs in dense and 
cold clouds Cg < 1 km s~^. Although these strong accretors emit a large X-ray luminosity, 
their Bondi - Hoyle radius ^ 10^^ cm can exceed the Stromgrcn radius. The effective 
accretion radius, where the gas is coupled through magnetic fields and angular momentum 
is removed by fast collisions, shrinks below tbh- Accreting NSs with rg < tbh may 
enlarge the cosmic ray (low FIP)/(high FIP) elemental ratio, that is observed to be ~ 3 
times larger than in Local Galactic matter (Meyer 1985). Selection between ions and 
neutral atoms occurs far from the NS, at the radius where each element becomes ionized. 
The radii of photoinization fronts would increase as the ionization potential decreases 
(note that the ionized continuum is much fiatter than blackbody). By assumption, ions 
are more effectively accreted than neutrals, and by such an ion - neutral bias the relative 
fraction of low FIP elements in the inflow is enhanced. Furthermore the ions would be 
anchored within the inflow by the frozen magnetic fields, while neutral atoms can diffuse 
in a direction oppose to the temperature gradient. As long as the diffusion time scale 
r'^/D is smaller then the dynamic time scale, r/v, ion-neutral separation by diffusion can 
be significant. The diffusion coefficient for neutrals in an ionized hydrogen gas (Geiss 
1982) i:> = 1.1 X 102iT/(10^ K)3/2n-^ cm^ s'^ may, indeed, be larger than vbhv ~ 
lO"*^® — 10"*^^ cm^ s~^. The ion-neutral bias is particularly applicable in dense clouds 77-4 ^ 1 
or, more precisely, when n\^^v^Q > 70, otherwise the accretion radius is well within the 
ionization sphere and the inflow is ionized in the flrst place. 

Gamma-ray observations show a sharp concentration of high energy radiation in a 
very thin Galactic disk, with evidence for imprints of the Galactic arms. If this emission 
comes from decaying pions, originating in cosmic ray interactions with the ambient gas 
(Stechcr & Stecker, 1970), it requires the cosmic ray concentration to follow the high energy 
photons. In our model this is fulflUed since the dense ISM clouds are also concentrated in 
the spiral arms, generally close to the galactic equator. 

What is the electromagnetic imprint of isolated NSs? The soft X-ray luminosity from 
the shock region is a fraction of oc Z/£.j.(l — fcr) J cind the expected spectrum is nonthermal. 
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From equations 13 and 14 one can predict a power law of spectral index (a — l)/2 ~ 
1/2. Other emission mechanisms, involved with heat loss of the star itself would basically 
produce thermal emission, but a magnetized atmosphere could affect the radiation transfer 
and modify the blackbody curve. A crucial point to note is that the Galactic diffuse X-ray 
background hardly restricts the cosmic - ray luminosity from isolated NSs. The integrated 
flux from isolated NSs was found to contribute no more than a fraction of a per cent to 
the soft X-ray background (BM93). The strongest accretors are located in dense clouds, 
and so even if they emit strong soft X-ray radiation, it would be greatly absorbed. Using 
a photoionization cross section ~ 10~^^(i?/KeV)~'^ cm^ (H — atom)~^, the cloud optical 
depth Tx ~ 10^(£^/KeV)~^n4(i?/10^°cm), where R is the cloud radius. Solitary isolated 
NSs, moreover, may be the origin of the enhanced X-ray emission within ±40*^ longitude 
of the Galactic center, as observed by EXOSAT (BM93, Maoz & Grindlay 1994). 

Although isolated NSs are the majority of the NS population, they have so far received 
only limited attention compared with that given to the brighter ones, the pulsars and the 
(NS) X-Ray sources. Only a small fraction of the Galactic isolated NSs can be observed 
directly. Following the early suggestion that nearby isolated NSs might be detectable due 
to radiation released from the accretion of the ISM (Ostriker, Rees & Silk 1970), Madau 
& Blase (1994) argue that the ~ 100 sources detected in the EUV - soft X-ray range the 
ROSAT - WPC and EUVE all - sky surveys, and which have no previous identification in 
the optical or other energies, are the natural candidates for isolated NSs. Recent obser- 
vations at the Wise Observatory of the fields of ~ 15 of these sources show that indeed 
in many cases the sources are likely hot or peculiar, and unlikely to be white dwarfs or 
late type fiare stars. Simultaneous observations planned from a number of experiments 
aboard the Spectrum X Gamma observatory, in UV {TAUVEX), EUV {EUVITA), and 
X-Ray bands {SODART, Jet-X, MART), will be extremely useful for the identification of 
such objects. The predicted characteristics to be studied in such UV - X-ray observations 
include spectral index, the presence of synchrotron absorption lines, temporal modulation, 
proper motion, and HII cometary - like regions. Detection of isolated NSs would be an im- 
portant step towards our understanding of this Galactic component. The model presented 
here can be constrained by observations of the total isolated NS population, their spatial 
and velocity distributions, and the fraction of them in dense clouds. 
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